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We show that azimuthal asymmetry in transverse momen- 
tum distribution of particles observed in non-central heavy 
ion collisions can originate from pre-equilibrium rescattering. 
The phenomenon is studied using computer simulation of ex- 
panding pion gas created in Pb+Pb 160 GeV/n collisions. 
Results obtained are discussed in comparison with experimen- 
tal results of NA49 collaboration. Theoretical understanding 
of mechanism generating the asymmetry is presented. Con- 
clusion about non-equilibrium properties of the phenomenon 
studied is drawn. 



I. INTRODUCTION 

Since the discovery of directed flow in pioneering ex- 
periments on Bevalac [[l] the phenomenon of azimuthal 
asymmetries has been studied extensively at various 
beam energies using different types of colliding nuclei. 
Besides the interesting phenomenon of directed flow two 
additional signatures of the collective behavior of hadrons 
have been identified in heavy ion collisions (H1C): 

1) Radial flow, existence of which has been clearly es- 
tablished in central collisions 0. Recently signature of 
radial flow in non- central collisions has been predicted 0] 
and most likely also confirmed experimentally 0j . 

2) Elliptic flow, which can be oriented orthogonally to 
the reaction plane (squeeze-out) at Bevalac energies pi or 
parallel to the reaction plane at AGS j(|[7j or SPS 
energies. In Ni+Ni 2GeV/n experiments orientation of 
the elliptic flow depends on transverse momentum region 
of neutral mesons selected for the analysis M. 

Phenomenon of azimuthal asymmetries is satisfactorily 
understood as a consequence of collective hydrodynami- 
cal behavior of nuclear matter described by the equation 
of state or as a result of shadowing effect and absorption 
processes in spectator parts of colliding nuclei. 

In this work we study properties of second-order el- 
liptic flow of final state pions. We concentrate mainly 
on anisotropy created in Pb+Pb 160 GeV/n non-central 
collisions studied recently by NA49 collaboration using 
Main TPC f§ and Ring Calorimeter || data. 

Second-order flow oriented parallel to the reaction 
plane was predicted in theoretical work of J.-Y.Ollitrault 
JlO| for ultra-relativistic energy domain. Mechanism 
of in-plane elliptic flow generation jlO) is based on the 
assumption of thermalization of nuclear matter in the 



asymmetrical overlapping region of non-centrally collid- 
ing nuclei. This assumption is essential for the existence 
of different pressure gradients in different directions in- 
side the asymmetrical overlapping region. During the 
subsequent hydrodynamical evolution different pressure 
gradients generate second-order asymmetry in transverse 
momentum distribution of particles ]ic| ]. Also in the re- 
cent theoretical work on elliptic flow at AGS energies 
an assumption about weak sensitivity of the result- 
ing flow asymmetry to non-equilibrium effects (in central 
rapidity region) is formulated. 

In this work we show that elliptic flow can be generated 
by rescattering process among produced hadrons (pi- 
ons) without the equilibration in prehadronic or hadronic 
stage of heavy ion collision. Our study is motivated by 
results of computer simulation p"l| of the rescattering 
process in the expanding pion gas created in Pb+Pb 160 
GeV/n non-central collisions. 

The paper is organized as follows: In section II we 
briefly describe our rescattering simulation of the ex- 
panding pion gas - the CASCUS 1.0 model. Main results 
of the simulation - centrality dependence, p t dependence, 
rapidity dependence and the equilibration dependence of 
the asymmetry are presented in section III. The results 
are compared with experimental data Q. Qualitative 
arguments offering our understanding of the mechanism 
generating the asymmetry are presented in Section IV. 
At the end of this work a summary of results is given 
and short conclusions are drawn. 



II. RESCATTERING PROCESS IN PION GAS 

Total number of pions generated in ultrarelativistic 
Pb+Pb collisions at SPS energy is surprisingly large. 
Even more than 2500 pions can be generated in a central 
Pb+Pb collision at this energy. Number of secondary 
particles thus exceeds number of primary nucleons in 
the system. Assuming that pions are created in inde- 
pendent nucleon-nuclcon collisions the initial transverse 
momentum and rapidity spectra of produced pions are 
determined mainly by the dynamics of pion production 
in nucleon-nucleon collisions []. 



*Even if this picture of pion formation may be oversimplified 
we believe that it can provide a reasonable approximation to 
reality for SPS energy region 
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Our rescattering simulation is based on the following 
scenario: Initial state of the pion gas is characterized by 
momentum and spatial distribution p(p, x, t) of pions as 
they are produced in nucleon-nucleon collisions. As the 
system evolves positions of pions change naturally due 
to the classical motion of pions. Momenta of pions can 
change as well in mutual collisions of pions. As the size of 
the cloud of pions increases, spatial density of pions be- 
comes smaller and consequently collisions among pions 
become rare. At some point collisions cease at all and 
momentum distribution of pions is not influenced any 
more. Pions move then freely on their way to detectors. 
Thus momentum distribution of pions recorded by detec- 
tors is not equal to the original momentum distribution 
of pions as they were produced in nucleon-nucleon inter- 
actions. It is changed during rescattering process in the 
early stage of the pion gas expansion. In next subsection 
we give short description of the computer simulation we 
have performed. 



A. Simulation of rescattering process 

Our rescattering model [jllj is similar to the rescatter- 
ing simulation described by T.J.Humanic p2[ . Two ad- 
ditional concepts have been used in our version of rescat- 
tering model: 

1) Initial state of the pion gas in non-central collisions 
was generated by cascade program |l3| and 

2) Formation time of hadrons H] was introduced. The 
simulation is performed in CMS of the Pb+Pb collision 
with time step AT = 0.1/m/c. In each time step po- 
sitions of pions are changed according to the classical 
equation of motion x(T + AT) = x(T) + v ■ AT. Two 
pions collide when their distance is smaller than critical 
distance d c (s) = y/ a(s)/it where a{s) is isospin averaged 
elastic pion-pion cross section (taken from experimental 
data flol ). New momenta of pions are determined in 
CMS frame of pion-pion collision using differential cross 
section of the elastic pion-pion interaction. Then new 
momenta of pions are transformed back to the global 
frame of the simulation (CMS of Pb+Pb system). 



After being created pions are not allowed to participate 
in mutual interactions for Lorentz dilated time interval - 
the formation time Tf = 7 • Tf and for interaction delay 
time Tj after each collision. Formation time Tf influences 
expansion of the pion gas mainly in the initial stage of 
time evolution, interaction time Tj suppresses unphysi- 



cal effects of rescattering simulations |16 and allows to 
influence total number of collisions in the pion gas. 

Since no collective behavior of primary nucleons is in- 
corporated in cascade generator (jl3| and since the pro- 
duction of pions in primary nucleon-nucleon collisions is 
azimuthally isotropic the initial state of the pion gas is 
azimuthally isotropic in momentum space also in the case 
of non-central collisions: 



g(x,p,t) = X (t)-A T (x)-S T (p) 



(1) 



Here A T denotes azimuthally asymmetrical and S T de- 
notes azimuthally symmetrical distribution in transverse 
plane. During the pion gas expansion symmetrical S T (p) 
distribution is changed by mutual collisions of pions. 

First results of our simulation (centrality dependence 
of the asymmetry) have been presented during Heavy Ion 
Workshop on Particle Physics (September 1996) in Slo- 
vakia [ pl]JT^Jl^ |. Detailed description of the simulation 
can be found in thesis |l9[] . In following subsections we 
present additional results of the rather extensive simula- 
tion we have performed. 

III. RESULTS OF THE SIMULATION 

Rescattering simulation was performed for random and 
unknown orientations of the reaction plane in events gen- 
erated by cascade program [ pl| |. During the subsequent 
analysis impact parameter orientation had been deter- 
mined for each event [fl7f and histograms of azimuthal 
distributions of pions had been filled (sec Fig. 2). Az- 
imuthal asymmetry of second order oriented in the re- 
action plane (not present in the initial state of the pion 
gas) had been clearly identified in transverse momentum 
distribution of pions after the rescattering process. 




Fig.l Time delay parameters incorporated in our rescat- 
tering simulation. 
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Fig. 2 Azimuthal distribution of pions before and after 
the rescattering process for events with b = 3fm and 7fm. 
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Strength of the asymmetry was evaluated by fit to the 
function: 

R(^) = S [1 + S 2 •cos(2<^)] (2) 

In next subsections we present main properties of the 
asymmetry found. 



We find also to be surprising, that the asymmetry is 
generated in the case of small number of collisions per 
pion (e.g 0.5 collision per pion). Dependence of the 
strength of asymmetry on the equilibration of the pion 
gas is studied at the end of this section. Now we present 
results on p t and rapidity dependence of the asymmetry. 



A. Centrality dependence 

Events with impact parameter 6 = 3, 5, 7, 9, 10, 11, 12fm 
generated by cascade program have been used as in- 
put for the rescattering program and subsequently an- 
alyzed for second-order asymmetry. Following behavior 
had been found for different values of rescattering model 
parameters r/,Tj: 



B. Transverse momentum dependence 

Transverse momentum dependence of the asymmetry 
has been identified in our results already in the early 
stage of this work jll] . It has been studied quantitatively 
for different values of rescattering program parameters 
Tf,Ti. 100 events with impact parameter 6 = 7fm have 
been used for this analysis. Selected results are presented 
in Fig. 4. 
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Fig. 3 Centrality dependence of the asymmetry for dif- 
ferent values of 17, T; parameters. Number of collisions 
per pion (indicated for each data point) decreases with 
increasing value of impact parameter. Non-zero asym- 
metry is generated also for small number of collision per 
pion (0.5 Coll/pion at b=llfm). 

Asymmetry increases with impact parameter up to the 
region 6 =9 -10 fm and then it decreases. This behavior 
has been predicted in the work of J.-Y. Ollitrault jl(| and 
it is in remarkable agreement with experimental data Q . 
From the point of view of rescattering simulation cen- 
trality dependence of the asymmetry can be understood 
easily. For 6 = initial conditions have cylindrical sym- 
metry and there is therefore no reason for the asymmetry 
in the final state. For very peripheral collisions interac- 
tions of pions are very rare due to low multiplicity of 
pions. Therefore asymmetry is expected to decrease for 
peripheral collisions. Sizeable asymmetry is generated 
in the intermediate region where spatial asymmetry in 
the initial state and number of pion-pion collisions are 
sufficiently large. 

Shape of the centrality dependence is not very sensitive 
to the number of collisions per pion in the expanding pion 
gas (see Fig. 3). 
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Fig. 4 p t dependence of the elliptic flow of pions for dif- 
ferent values of parameters Tf,Ti. 

The asymmetry increases with p t for number of col- 
lisions per pion N co u / w = 4, 8, 14 in rough agreement 



with experimental data W. For N t 



coll. /it 



0.6,0.3 the 



shape of the p t dependence changes together with abso- 
lute strength of the asymmetry. Decrease of the asym- 
metry strength at large p t for N co n/ n = 0.3,0.6 can be 
understood easily. High p t pions can escape from the 
pion gas volume without interaction due to large value of 
Lorentz-dilated formation time Tf and also due to their 
larger velocity. Therefore smaller asymmetry is gener- 
ated for rarely colliding high-j^ pions while larger asym- 
metry is generated for colliding low-pt pions. Decrease 
of the asymmetry in case N co u j v = 8, 14 in comparison 
with N co u_/ n — 4 is discussed in section D. 

Detailed comparison of pt dependence of asymmetry 
with experimental data may allow to determine physical 
parameters incorporated in the simulation (in particular 
the formation time). 

For this goal however improvements in our rescatter- 
ing model (e.g. dynamics of Ti delay parameter) and 
also higher statistics of experimental data Q would be 
necessary. 
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C. Rapidity dependence 

Rapidity dependence of the asymmetry has been stud- 
ied using b = 7fm events sample. Resulting rapidity 
dependence found by our simulation exhibits maximum 
at distance Ay = 0.75 from central rapidity region. We 
find to be interesting that results on rapidity dependence 
of elliptic flow presented during Quark Matter 97 confer- 
ence j|] also exhibit decrease of the asymmetry strength 
at mid-rapidity region. 
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Fig. 5 Rapidity dependence of the asymmetry. Note de- 
crease of the asymmetry at midrapidity. 

This result is worth of further investigation. It is not 
excluded that the origin of the mid-rapidity decrease of 
the asymmetry is related to behavior studied in the next 
subsection. 



D. Thermalization dependence 

Numerical value of the rescattering model parameters 
(see Fig.l) allows to influence total number of collisions 
in the expanding pion gas. We have studied dependence 
of the asymmetry on total number of collisions in the 
pion gas for the set of b = 7fm events. Following result 
has been obtained (Fig. 6): 

With increasing number of collisions in the pion gas the 
asymmetry first increases, it exhibits maximum close to 
Ncoii/x = 4 region and then it decreases. This decrease 
depends mainly on number of collisions per pion, it is 
not a consequence of the physics of formation time Tt or 
interaction time Tj parameter. Both Tf = 0.5, Tj = 0.2 
and Tf = 0.2, Ti = 0.5 combinations lead to decrease of 
asymmetry from Tf = 0.5, 1$ — 0.5 level (see Fig. 6). This 
result can be interpreted as a non-equilibrium feature of 
the mechanism generating the asymmetry. We think that 
monotonic increase and saturation should be expected if 
the asymmetry were generated in an equilibrium process. 
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Fig. 6 Thermalization dependence of the asymmetry. 
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asymmetry seems to saturate. We guess this corresponds 
to hydrodynamical limit studied in work JTo| . 

At this point we can try to understand the decrease 
of asymmetry at mid-rapidity found in our simulation 
(Fig. 5) and present also in experimental data The 
decrease can be explained as a consequence of smaller 
number of collisions per pion (N co uj^ < 4) at midrapid- 
ity in comparison with number of collisions per pion at 
the maximum of rapidity dependence^} This point re- 
quires further and more careful study since the position 
of the maximum is different in our results and experi- 
mental data 0. 



E. Time evolution of the asymmetry 

Numerical simulation of the pion gas expansion allows 
us to study also quantities which cannot be accessed from 
final stage of the pion gas. We have studied time evolu- 
tion of the number of collisions in the pion gas together 
with the time evolution of second-order asymmetry. For 
the sake of simplicity the asymmetry has been character- 
ized by R p parameter: 



P (Pi) 



(3) 



which does not require the fit procedure for its evaluation. 

From results shown in Fig. 7 we conclude that the 
asymmetry is generated during the early stage of the ex- 
pansion. This confirms our conjecture that the asym- 
metry can originate during pre-equilibrium stage of the 
pion gas expansion. Theoretical understanding of mech- 
anism the asymmetry is generated by is presented in next 
section. 



2 This explanation is supported by the simulation of dilepton 
production from pion gas [EOl. 
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Fig. 7 Time evolution of the symmetry parameter R P (t) 
and evolution of number of collisions during the expan- 
sion of pion gas. Time dependence of R P (T) is nearly 
identical with the behavior of Sa(T) parameter [JL9| . 



Asymmetry of collision integral in Eq.(^) originates 
from the asymmetrical spatial distribution of pions in 
transverse plane. We do not attempt to evaluate the col- 
lision integral here though we guess this is accomplish- 
able task. The asymmetry of the collision integral can 
be understood intuitively using Toy model of pion gas 
expansion described in the previous paper pjj . In next 
section we derive formalism of two point decomposition 
of source which (we think) allows to study the asymmetry 
theoretically. 

A. Two point source of pions. 



IV. THEORETICAL ANALYSIS OF 
ASYMMETRY 

Expansion of the interacting pion gas can be described 
by kinetic equation which (in the case of zero external 
macroscopic field) takes the form: 



df(p,x, t) 

at 



V ? f{p,x,t) = [St 



Coll 



(4) 



\St\ Coll is the collision (or scattering) term describing 
influence of pion-pion collisions on distribution function 
f(p,x, t). 

In slightly simplified case (we do not take into account 
time interval during which pions are produced) our initial 
state of the pion gas - Eq. (Ty) can be written in the form: 



p(x,p,T ) = A 7 \x u x z ) ■ S 7 \ptiPz) 



(5) 



where again A T (x t , x z ) denotes asymmetrical spatial dis- 
tribution of pions in transverse plane and S T (pt,p z ) is 
symmetrical distribution of pions in transverse momen- 
tum. During the process of pion gas expansion asymme- 
try in spatial distribution of pions "leaks" into momen- 
tum distribution of pions which becomes asymmetrical 
in transverse momentum space. 

This phenomenon results from mutual interactions of 
pions in the expanding asymmetrical volume. If there 
were no collisions of pions, momentum distribution of 
pion gas (momenta of pions) would be unchanged. This 
is clear intuitively and it can be seen also on analytical 
level after the integration of kinetic equation ([|) in spatial 
coordinates and after the subsequent time integration: 



Fa(PuPz,T) = Fs(p t ,p z ,T ) 



[St] d 3 xdt (6) 



Here distribution F(p,t) = J F(p,x,t)d 3 x is distribu- 
tion of pions in momentum space. Since original momen- 
tum distribution of pions Fs(p,Tq) is symmetrical in p t 
plane (our initial condition) and since resulting distribu- 
tion Fa (p, T) is asymmetrical (this distribution is mea- 
sured by detectors) the collision integral f T J [St] d^dt 
must be asymmetrical. 



Let us consider two point-like sources of pions sepa- 
rated by distance D in transverse plane. Each source ra- 
diates pions with random azimuthal angle and (for sim- 
plicity) with fixed absolute value of momentum. If we 
switch off rescattering of pions emitted from these two 
sources azimuthal distribution of pions observed in a dis- 
tant detector surrounding the source (D <C R) is con- 
stant - symmetrical. However if collisions of pions are 
allowed the resulting distribution of pions is asymmetri- 
cal in transverse momentum plane. Excess of pions in 
direction orthogonal to D is generated as a consequence 
of collisions (see Fig. 8). 





Fig. 8 Azimuthal asymmetry generated by interactions of 
pions emitted from two point-like sources of pions. 

Elliptic asymmetry can be characterized by R p param- 
eter - Eq.(ra) and therefore we shall express pf, + p 2 , 
quantity for a given pair of pions after the collision at 
any place of transverse plane. If p 2 , + p 2 , has increased 
then also R p parameter increases what means that ellip- 
tic asymmetry is generated. After Lorentz boosting back 
from CMS of 7T7t collision to the rest frame of two-point 
source and after averaging over random orientations of 
pion momenta after the collision (in CMS of tttt collision) 
we obtain the result (see Appendix): 



(pi 



Pu 



Pcms + 7 2 /3y ■ MPCMS, P) 



2E CMS - 



Pcms 
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7+1 



(7+T? 



(7) 



and 



where A(p C MS,0) = 
(3 = J£q^g| ■ For {p 2 x , +p 2 ,,) a the same expression can be 
found just (3 y is replaced by (3 X . Increase of R p parameter 
is generated by positive values of 
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Ap' = (p: 



- P 



in 



(8) 



quantity (p yi 

i Py2 

momenta of pions before collision) 
and negative values of similarly defined Ap 2 quantity. 
Numerical values of Ap 2 and Ap 2 depend on coordinates 

of the collision point in transverse plane since (3,pcMS 
and also pi,P2 depend on coordinates of the collision 
point. In Fig. 9 we show Mathematica |22| contour-plot 
of Ap 2 quantity for two point source radiating pions with 
equal (pi = P2) and different (pi > p-i) momenta. 

Two-point source represents (in rough approximation) 
overlapping region of colliding nuclei in transverse plane. 
It is clear, that resulting elliptic asymmetry asymmetry 
is oriented in the direction of impact parameter - orthog- 
onally to D vector. 





x[fm] 



x[fm] 



Fig. 9 Contour plot of Ap y - see Eq.(^) for a two point 
source emitting pions with momenta pi and p2- White 
color represents positive values of Ap y quantity. 

In next section we study asymmetry generated by non- 
trivial discrete and also continuous source. 



B. Two-point source decomposition 

Main idea of the two-point decomposition of source is 
following: Each spatially distributed source of interacting 
particles behaves (in 1 collision/pion approximation) as 
a set of two-point sources (see Fig. 10). Each two-point 
source (characterized by vector D) generates asymme- 
try with the orientation determined by vector D. Our 
intention is to express analytically resulting asymmetry 
generated by the whole source. 

Let R2{4>, D) is azimuthal distribution of pions (parti- 
cles) emitted from two-point source characterized by vec- 
tor D. Asymmetry of R2(4>, D) is generated by collisions 
of pions at any place of transverse plane. Contribution 
from all possible collision points is averaged. Let us also 
have non-trivial e.g. discrete source radiating indepen- 
dent pairs of interacting pions. (This is our specification 
of 1 collision per pion approximation.) 
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Fig. 10 Discrete asymmetrical source of pions and its 
two-point decomposition. Note the excess of horizontally 
oriented two-point sources in comparison with vertical 
two-point sources. 

Then in lcoll./7r approximation resulting azimuthal 
asymmetry of particles emitted from more complex but 
discrete source (see Fig. 10) is given by weighted sum of 
azimuthal distributions of all possible two-point sources: 



(9) 



Weights Wij depend on intensity of point-like sources 
i, j. Expression (||) can be extended to the case of spa- 
tially continuous source. 

If we define probability distribution for the emission of 
a pair of particles from any places X\,X2 in the source 
p(x) at relative distance D = x\ — X2 



W P (D) 



d 2 x\d 2 X2 



p(xi)p(x 2 ) ■ 8{D ~ (xi - x 2 )) 

(10) 

what can be rewritten as 

,x + D 



WJD) 



(11) 



then resulting azimuthal distribution of particles emitted 
by spatially distributed source p{x) can be expressed in 
1 collision per pion approximation as follows: 



Rptt) = / W p {D)-R 2 {^D)d 2 D 



(12) 



Dependence of the resulting R p ((j)) distribution on e.g. 
formation time and 7T7t cross section is encoded in distri- 
bution R2(cj>,D). If properties of i?2(</>,-D) distribution 
are known, expression @ allows us to say which source 
can generate asymmetry at analytical level - without nu- 
merical simulation. 

Asymmetry of resulting R p distribution is generated 
only by asymmetrical part of W p (D) distribution 



A(D,<fi) = W p (D,<p)-S(D) 



(13) 



where S(D) is symmetrical part of W p (D,(j>). S(D) is 
equal to minimum of W P (D, (f>) at given D. Asymmetri- 
cal source distribution p(x) leads to asymmetrical W P (D) 
distribution which results in non-zero asymmetrical dis- 
tribution A(D,(f)). Non-zero A(D,(p) distribution gener- 
ates asymmetrical R p {<j)). Properties of R p (4>) depend on 
analytical properties of A(D,(j>) and R2(4>,D). 
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V. SUMMARY AND CONCLUSIONS 

We have shown that second order azimuthal asymme- 
try of hadrons known also as in-plane elliptic flow can be 
generated by rescattering process. We have performed 
extensive study of this asymmetry using computer simu- 
lation of Pb+Pb 160 GeV/n non-central collisions. Main 
properties of the asymmetry (impact parameter depen- 
dence, pt dependence and rapidity dependence) studied 
by the simulation are in rough agreement with experi- 
mental data 0] . 

Transverse momentum dependence of the asymmetry 
(Fig. 4) is found to be sensitive to formation time param- 
eter. Careful analysis of experimental data together with 
improved simulation of hadron gas expansion might allow 
to determine experimental value of the formation time of 
pions from data. 

Rapidity dependence of the second-order asymmetry 
of pions (Fig. 5) exhibits minimum at midrapidity. Study 
of thermalization dependence of the asymmetry (Fig. 6) 
allowed us to interpret this behavior as a consequence 
of different number of collisions among pions in different 
rapidity intervals. 

Dependence of the strength of asymmetry on total 
number of collisions in the pion gas (Fig. 6) together with 
non-zero value of resulting asymmetry in the N co u/ 7T =0.5 
case (Fig. 3) instigates us to conclude that the studied 
asymmetry can be generated in a non-equilibrium (pre- 
equilibrium) process. This conclusion substantially mod- 
ifies our point of view on physical meaning of the in-plane 
elliptic flow of pions measured in experimental data. 

We have developed formalism of two-point source de- 
composition (in one collision per pion approximation) 
which allows theoretical study of the asymmetry. We 
have shown that resulting distribution of pions in mo- 
mentum space is related to the initial spatial distribution 
of pions. 

We conclude that azimuthal anisotropy in momen- 
tum distribution of hadrons can be generated by pre- 
equiUbrium rescattering process. 



VI. APPENDIX 

For a given pair of pions colliding at place (x, y) of 
transverse plane the absolute value of momenta in CMS 
frame of collision is: 



Pcms = \/M(x,2/) 2 /4 - ml 



(14) 



where M(x, y) 2 = {E 1 +E 2 ) 2 - (pi(x, y) +p 2 (x, y)f . Let 
the orientation of momentum p CM g after the collision in 
CMS of 7T7T collision is random and characterized by angle 
a: p^ MS = p CMS ■ cos(a);p^ MS = p CMS ■ sin(a). Then 
momentum of pion after the collision in the LAB frame 
(rest frame of two-point source) is: 



vi = ?; cms 



+ 7/3'(-^ T /3'-P 1 
7+1 

V1+P2 



'CMS 



E, 



CMS! 



(15) 



where (3' = —f3 = — jg*^ • Consequently for e.g. y 
components of pion momenta after collision we have: 



Pyi 
Py 2 



Py^+lPyG, 



+ iPy ■ G 2 



'CMS 



(16) 



where G\ 



CMS 



Ecms an d similarly G 2 



,'2 



Pyi +Pm 



Ecms- Then the sum p y 2 ± + p y ~ gives: 

2p 'CMS 2+7 2p2 Gt 2 +2p >CMS 10yGl2 



where Gj~ 2 = G\—G 2l G{ 2 = Gf + G 2 and relationship 
^cms _ _p^CMS been used. Equation (^) is valid 
for any given angle a of pion momentum after the colli- 
sion in CMS. For the sake of simplicity we assume that 
this angle is random what allows to perform averaging of 
result ([l7]) over values of a. This leads to the result: 



/ 2 1 2\ 2 1 2 o2 

\Pyi +Py2)*= Pcms +7 P y 
27 2 



+ 



7 + 1 



PCMsfi 



{^-^) 2 P 2 PCMS+ 2E CMS 
(18) 



where following sub-results have been used: 

7 2 /3 2 



{G{ + G 2 2 ) a = 2E, 



CMS' 



(7 + 1) 

(sin a ■ [Gi - G 2 }) a = 2p C MS- 



2 Pcms 



7 



7 + 1 



(19) 
(20) 



Result (|T^) can be rewritten directly in the form of 
Eq.@. 
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